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Abstract 

Following our recent letter [T], we study in detail an entry-wise diffusion of non-hermitian complex 
matrices. We obtain an exact partial differential equation (valid for any matrix size N and arbitrary 
initial conditions) for evolution of the averaged extended characteristic polynomial. The logarithm of this 
polynomial has an interpretation of a potential which generates a Burgers dynamics in quaternionic space. 
The dynamics of the ensemble in the large N is completely determined by the coevolution of the spectral 
density and a certain eigenvector correlation function. This coevolution is best visible in an electrostatic 
potential of a quaternionic argument built of two complex variables, the first of which governs standard 
spectral properties while the second unravels the hidden dynamics of eigenvector correlation function. We 
obtain general large N formulas for both spectral density and 1-point eigenvector correlation function valid 
for any initial conditions. We exemplify our studies by solving three examples, and we verify the analytic 
form of our solutions with numerical simulations. 

Keywords: Non-hermitian random matrix model. Diffusion equation. Burgers equation. Characteristic 
polynomial 

PACS: 05.10.-a, 02.10.Yn, 02.50.Ey, 47.40.Nm 


1. Introduction 

The concept of matrices filled with entries subject to the diffusion process was first introduced by 
Dyson [2] and applied in the context of both Gaussian Unitary Ensemble (GUE) and Circular Unitary 
Ensemble (CUE). The arising Coulomb gas analogy had a major impact on understanding of random ma¬ 
trices [3]. Today the marriage of stochastic processes and random matrices still brings several new insights. 
Examples include the study of determinantal processes msiis], Loewner diffusion [7] or the fluctuations of 
non-intersecting interfaces in thermal equilibrium [5]. 

Recently, several authors B iin] have approached the diffusion in the GUE from a new perspective. 
They found a viscid complex Burgers equation for the logarithmic derivative of the averaged characteristic 
polynomial fjg{z,T) = ;^i9z In {7Ar(z, t) (associated with a Hermitian matrix filled with entries performing 
Brownian motion in the complex space): 

drfN{z,T) + fN{z,T)djN{z,T) = “(z, t), (1) 
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where r is the diffusion time and z is a complex variable. The role of viscosity is played by the inverse 
size of the matrix N. In the —>■ oo limit, becomes the Green’s function G{z,t) and the partial 

differential equation becomes inviscid: 

drG{z,T) + G{z,T)dzG{z,T) = 0. ( 2 ) 

A solution of the latter equation (by the method of characteristics) requires an introduction of shocks, which 
turn out to coincide with the edges of the spectra. This phenomenon leads to a novel interpretation of known 
matrix results, since microscopic universal behavior of the spectra emerges as an expansion around the shock 
wave of the viscous equation. Nontrivial initial conditions give rise to shock collisions which are equivalent to 
the merging of the spectrum boundaries. This way, not only Airy but also the Pearcey functions are captured 
in the same formalism. A similar Burgers equation was also obtained for the Wishart ensemble and chiral 
GUE yielding a universal scaling associated with the Bessoid function m- The equivalent phenomenon 
appears also at the level of CUE diffusion, providing new insight for order-disorder transition of Wilson 
loops in Yang-Mills theory [12]. 

Recently, this program got extended to the realm of matrices with complex eigenvalues [1]. The success 
of such an extension is a priori surprising, since hermitian and non-hermitian matrix models seem to be 
hardly comparable. In the former, the hermiticity condition confines the eigenvalues to the real axis. In 
the latter, there are no constraints, and the eigenvalues spread over the whole complex plane. Furthermore, 
non-hermitian models develop discontinuities at the spectral density boundaries, a feature observed even for 
the well-known Ginibre Ensemble [13], for which the spectral density is given by: 

p{z,z) = \l-\z\\ (3) 

TT 

whereas e.g. in the GUE case, the Wigner semicircle is continuous across the spectral edge, and only the 
derivatives are discontinuous. Other differences arise when one considers the evolution of the diffusing 
matrices. In the hermitian case, the evolution is determined by the initial eigenvalues only, whereas in the 
non-hermitian models, the information on initial eigenvectors additionally affects the shape of the spectral 
density. 

This paper is a continuation and an extension of the ideas of non-hermitian diffusion announced briefly 

in HI- 

Our approach is rooted in the standard electrostatic analogy, but requires a novel setting (we call it 
the quaternionic method), which can be viewed as an extension of the standard Dysonian strategy applied 
originally to the hermitian (or unitary) matrix diffusion case. The main objects of our interest are the 
spectral density (obtained from the resolvent) and a certain one-point eigenvector correlation function. We 
stress that the aforementioned eigenvector correlator is crucial for understanding the diffusion process of non- 
hermitian matrices. We also point out, why the importance of this correlator was disregarded in majority 
of the studies of non-hermitian random matrix models. 

The basic object of our studies is an averaged “extended” characteristic polynomial (AEGP). An exten¬ 
sion follows from an introduction of two pairs of complex variables (compared to one complex variable in 
standard treatments). Surprisingly, AECP obeys a certain partial differential equation akin to the diffusion 
equation, for arbitrary size of the matrix and for arbitrary initial conditions, and is exactly integrable. The 
diffusion happens in the auxiliary plane “perpendicular” to the complex plane where the eigenvalues reside. 
In the large matrix size limit, the logarithm of the AECP can be viewed as an electrostatic potential and its 
derivatives with respect to the two complex variables yield a pair of coupled Burgers-like equations for the 
non-hermitian Green’s function and eigenvector correlation function. We would like to mention, that in the 
standard electrostatic analogy [HUSIIIS] the ’’second variable” is treated as an infinitesimal regulator only. 
This is the reason why the dynamics, as a function of this variable, remained hidden, and the complementary 
information on the eigenvector correlator co-evolving with the spectra was absent. 

To illustrate our findings, we consider a couple of examples of initial conditions. In most of them 
we demonstrate the explicit solutions of Burgers equations, obtaining formulas for the spectral density, 
eigenvector correlators and the electrostatic potential. We note that an inspection of the Burgers-like 
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equation identifies the shock line with the non-holomorphic sector of the spectral density. Moreover, we 
show the insensitivity of the shock formation to the initial condition chosen. This hints to a lack of truly 
new universality classes in this type of models, which is furthermore corroborated in the study of the universal 
behavior in the vicinity of the spectral collision in one of the examples. 

The paper is organized as follows. In Sectionwe define Dysonian non-hermitian diffusion. In Section]^ 
we briefly review the electrostatic analogy and the quaternionic method. We proceed in Sectionj^by deriving 
the partial differential equation for the AECP and presenting its integral representation. In Section we 
derive a pair of coupled Burgers equations for the on- and off-diagonal parts of the quaternionic Green’s 
function (in the large N limit) thus making a link to the quaternionic method. Subsequently we solve 
them with the method of complex characteristics. Finally, we obtain an implicit solution to the equation 
for the potential in terms of the Hopf-Lax formula and derive large N formulas for the spectral density, 
eigenvector correlation functions and the boundary of the spectrum valid for an arbitrary initial matrix. 
Section is devoted to the examples of a) Ginibre, b) spiric and c) I-band non-normality matrices. We 
apply previously described methods to these cases, depict the characteristics picture and obtain the large 
N limit spectral density and eigenvector correlators. We also comment on critical behavior of the EAGP. 
Section is devoted to a curious observation by Osada which actually has triggered our interest in the 
diffusion of the Ginibre ensemble. We provide also an explanation of the Osada’s observation. Section 
summarizes the paper and outlines some possibilities of further investigations of the observed patterns. 

Three appendices hide technicalities: Appendix A demonstrates the derivation of the key diffusion equa¬ 
tion for an AEGP, Appendix B clarifies the link between Ginibre and Wishart/chiral ensembles, and Ap¬ 
pendix G determines weights and normalizations needed for establishing the universal scaling at the shock. 

2. Dysonian non-hermitian diffusion 

Gonsider a non-hermitian N x N matrix X = whose elements Xij = Xij -I- iyij 

undergo 2A^^-dimensional Brownian motion 

(r) = (r), (r) = dB^ (r), (4) 

where and B^^ are independent Wiener processes. We restrict ourselves to deterministic initial conditions, 
that is we assume that each element of the matrix has a given initial value Xij = ixo)ij and yij = {yo)ij for 
T = 0. This condition can be concisely written as X = Afo, for t = 0. Glearly, this model is a straightforward 
extension of the Dyson random walk [2] to the realm of non-hermitian random matrices. 

The joint probability density for matrix elements evolves according to the 2A^^-dimensional diffusion 
equation 


drP{X,r) = ^ + 91,)P{X,t), (5) 

ij 

with the initial condition P{X, 0) = 5{X — AIq). The probability measure for random matrices at time r is 
defined by dy,r{X) = 'D[X]P{X, r) where T’[A'] = dxijdyij^ and the statistical averages by 

{F(X)), = J V[X]P{X,t)F{X). (6) 

The hermitian version of the model, discussed by Dyson, reduces to a model of evolution of eigenvalues. In 
that case eigenvectors can be integrated out. What makes the non-hermitian extension interesting is that 
in addition to eigenvalues one has to control also the evolution of eigenvectors. We present a systematic 
method to do so. 
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3. Electrostatic analogy and quaternions 


In this section we briefly recall the method to calculate eigenvalue distribution of random matrices in 
the limit N ^ oo. The method is based on “electrostatic” analogy [IHIIlllle]. One defines a quantity 

^{z,w,t) = ^(Trlog((z-X)(z-X^) + |wp))^, (7) 

which can be interpreted in the limit ic —)■ 0 as an electrostatic potential of a cloud of N electric charges 
interacting on the ^-complex plane. The corresponding electric field is 


G {z, w, t) = dz^iz, w, t) 


1 

iv 



z-Xt 

{z-x){z-xr) + \w\^ 


( 8 ) 


Identifying the real and imaginary part of as vector components G = {E^ — iEy)/2 one can rewrite the last 
equation in the vector notation as E = (Ex, Ey) = Vz®. The minus sign in front of Ey and the scale factor 
1/2 in the relation of G to electric field if is a matter of convention. 

We are interested in the eigenvalue distribution 


= (z-Zi\ , 


(9) 


where ZiS are the eigenvalues of X. The limiting eigenvalue density can be calculated from the Gauss law 


p(z, t) = —dzG (z, w,t) , w — 7> 0. 


( 10 ) 


This relation follows from a standard representation of the complex Dirac delta function (z — Zi) = 
r. The expression in the brackets on the r.h.s. of (81 can be cast into the standard form 


Nr 


lim ■- 

Itul-s-o (Nr+b-^il ) 
of resolvent (z — at the price of introducing 2N x 2N matrices 


Q = 


z —w 

w z 


x = 


X 0 

0 xf 


in place of the original N x N ones. The resolvent is a 2 x 2 matrix 

where the block-trace is defined as 


( 11 ) 


( 12 ) 


bTr 




fTiA 

TrB\ 

[g 

D 

~ iTrC 

Tr£»j 


We refer to Q(z,w,t) (12) as to generalized Green’s function or quaternionic resolvent [HI [19]. Note, that 


we use the representation of the quaternion in terms of Pauli matrices, i.e. Q = qol 2 + * 

-z = Qo + iqi and —w = q 2 +iq 3 - The diagonal element of the quaternionic resolvent Qn is equal to G(z, w, r) 
([^. The extended Green’s function Q{z,w) is an advantageous object since one can apply geometric series 
expansion to (Q — X)~^ which has a diagrammatic interpretation. This leads to a closed set of Dyson- 
Schwinger equations enumerating underlying planar Feynman diagrams. From these equations one can 
derive an exact form of the Green’s function in the limit X —>■ oo, as well as matrix-valued addition and 
multiplication laws |20j . We mention that the quaternionic extension is equivalent to another approach 
known under the name of hermitization method [SHEllEl!, in which the diagonal and off-diagonal blocks 
of matrices Q and X are flipped before the block-trace operation. 
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Having determined the quaternionic resolvent Q(z,w,t) one can determine the potential ^{z,w,t) or 
vice versa, since the two objects are related by a simple relation: 


G = 




(13) 


As follows from (10), the eigenvalue density can be derived from the potential ^{z,w,t) using the Poisson 
equation 


p{z,t) = -(92z$(z,0,r). 

TT 


(14) 


It turns out that the potential $(z,w,r) encodes also information about the correlations of eigenvectors 
|241 . being the special case of the Bell-Steinberger matrix naiiaisT]- One defines the correlation function 


0{z,t) = 




('^ ^a) 


(15) 


with Oap = {La\Lp) {Ra\Rp) where |Lq) (|-Rq)) are the left (right) eigenvectors of matrix X. It can be 
shown [28] that, in the N ^ oo limit, this correlation function is related to the off-diagonal elements of the 
resolvent as 0{z,t) t)Gii{z,0,t). Applying (^13^ we have 


0(z,t) = - d^<^(z,w,T)dtu<^{z,w,T)\^^^ = -|y(0,O,r)|2, 

TT TT 


(16) 


where V(z,w,t) = 9u,4)(z, w, r) is the velocity field, which plays the same role in the rc-complex plane as 
the electric field G{z,w,t) = dz^{z,w,T) in the z-complex plane. It is a vector field. If we parametrize 
positions on the w-complex plane as w = a + ib and V = (14 — iVb)j2, then V = Vu,$. The term ’’velocity” 
is related to the underlying Burgers dynamics to be discussed later. To summarize, the limiting eigenvalue 
density and the eigenvector correlation function can be calculated from the electrostatic potential using eqs. 


(14) and (16), respectively, and taking the limit w —>■ 0 which project quaternions to the z-plane. 

It remains to show how to calculate the electrostatic potential $(z,w,r) ([^. The standard method 
is based on enumeration of planar diagrams as mentioned above. The main object in this method is the 
Green’s function G- 

Here we propose an alternative approach which is based on the diffusion equation in the quaternionic 
(2 -|- 2)-dimensional space in the direction perpendicular to the z-complex plane. The primary object in this 
calculation is an extended characteristic polynomial defined in the next section. 


4. Averaged extended characteristic polynomial 

In order to calculate the potential ([^ in the limit A —>■ oo we rewrite it as 


1 


$(z,w,t) = — (logdet ((z - X){z - X'^) + . 


N 

and define an associated object - an effective potential 
(j){z, w,t) = ^ log (det ((z - A)(z - . 

where 


D{z, w, t) = (det ((z — X){z — X^) + = ( det 



(17) 

= ^logn(z,ri;,T), 

(18) 

fz-X -w \\ 

1 ^ -^-xy/r' 

(19) 


^Note that we introduced an additional 1/A^ factor as compared to the definition given in m to obtain a limiting density. 
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We refer to D{z,w,t) as to averaged extended characteristic polynomial (AECP). For the Gaussian process 
Q the determinant D{z,w,t) self-averages for TV —>■ oo and $ can be replaced by (j) in this limit. The 
advantage of using the latter is that the averaged extended characteristic polynomial D{z,w,t) (19), which 
appears in the definition of ^ ^ log D obeys a simple diffusion equation with respect to the variable w 


drD{z,w,T) = —dn,wDiz,w,T), 


( 20 ) 


as shown in Appendix A. Note that, from the point of view of this equation, z is a dummy parameter. The 
z-dependence appears only in the initial condition 


D{z, w, 0) = Do{z, w) = det (^{z - Xo){z - -b , 


( 21 ) 


which is completely determined by the initial matrix Xq. In other words, for each z we have an independent 
diffusion in the perpendicular rc-complex plane. The problem is therefore exactly integrable, and the solution 
of the diffusion equation (20) reads 


D(z,w,t) = ^ J exp^-A^ *"^ Do{z,w')d^w' , 


where Do(z,w') = D(z,w',0). The solution can be equivalently written as 


1 N f ( 

())(z,w,r) = — log— / ex^pN [cl)o{z,w') - 

N TTT Jc \ 


\w — w 


/|2 


d^w\ 


( 22 ) 


(23) 


where (j)o{z,w') = ^ logDo{z,w'). In the limit N ^ oo the last equation assumes the form of the Hopf-Lax 
formula 1^ 


(f){z, w, t) = max ( (j)o{z, w') — 

w' y 


\W — W 


(24) 


This equation describes the evolution of the electrostatic potential in the limit —)■ oo for the given initial 

configuration (f>Q{z,w), so this equation solves our original problem. 

A few remarks are in order. The characteristic polynomial D{z^w,t) and the potential (j)(z,w^T) (18) 
depend on w only through the norm |wp. The diffusion preserves the spherical symmetry of these quantities 
in the w-complex plane, so it is convenient to rewrite these equations in the radial part r oiw = re®“ skipping 

Do{z,r')dr', 


the dependence on the phase a. In particular (22) takes the form 

^2 _|_ j,/2 


r.. , 2N , f (2NTr'\ 

Diz,r,T) = —J^ r exp(^-iV—— 


and (24) simplifies to 


^(z, r, r) = m^ax (^q{z, r') - 


(r - r'f 


(25) 


(26) 


The second remark is on the role of the parameter w. Originally it was introduced as a regulator to the 
expression for the potential ([^ and eventually sent to zero. Here we promote w to a full complex variable 
and analyse the dynamics of the model on the entire rc-complex plane. This approach allows one to trace 
not only eigenvalues but also eigenvectors of the random matrix X and to break the symmetry between 
matrices having identical eigenvalues but different eigenvectors. 

A complex valued matrix can be Schur decomposed X = U{A + T)U^ where 17 is a unitary matrix, 
A is a diagonal matrix containing the complex eigenvalues, and T is a strictly upper-triangular matrix 
encoding information about eigenvectors. Two different matrices Xi,X 2 having the same eigenvalues A but 
different eigenvectors have different Ti and T 2 . The averaged extended characteristic polynomial (19) for 
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these matrices differs Di(z,w) ^ D 2 {z,w) when |ui| 7 ^ 0. From this difference one can read off information 
about eigenvectors. 

Moreover, the dynamics of the model in the quaternionic space has a beautiful physical interpretation 
in terms of the Burgers dynamics. The behavior of the model on the z-complex plane is a shadow of this 
dynamics. In particular, the support of the eigenvalue density p{z, t) coincides with the location of shocks 
of the quaternionic Burgers dynamics in the full quaternionic (z, r/;)-space. 


5. Burgers dynamics 


In the previous section we have found a solution to the diffusion equation (20). Here for completeness 


we relate the diffusion equation to Burgers dynamics. Using the definition (18) it is easy to see that the 


effective potential (j) = (j){z,w,T) and its gradient v{z,w,t) = dy^(j){z,w,T) fulhll the following differential 
equations 




(27) 


and 


^ ^WW T 1^1 ; 


(28) 


respectively. These equations describe Burgers dynamics on the w-complex plane for a two dimensional 
velocity field v = (14 — *I4)/2 derived from the potential </>: U = = {da4>, db(j)) where w = a + ib. The 

coefhcients of the Laplacian term can be identified as a hydrodynamic viscosity parameter = ^. Equation 
(28) follows from (l20| by an inverse Cole-Hopf transformation [30]. One can also write an equation for the 


z-gradient, g{z^w,T) = dz4>{z,w,T)\ 


dr9 = + d^\v\^. (29) 

The two gradients are related to each other as d^g = d^v. The effective potential (j) and the gradients 
reproduce the electrostatic potential and the quaternionic Green’s function in the limit IV —)■ 00 




g V 
-V g 


g. 


Let us now discuss the inviscid limit N —>■ 00 . The effective potential (27) obeys the equation 

which after applying rotational symmetry (in the variable w) simplifies to an equation 

dr(l) = ^(3r<(>)^ 


(30) 


(31) 


(32) 


for the radial variable r = |w|. The solution is given by the Hopf-Lax formula (26) which in our case is 
equivalent to 


(()(z,r,r) = (j)o{z,r^) - 


{r - Y 


(33) 


with r* being the location r* = r' of the maximum (26) given by the usual extremum condition 

o ^ ^ \ 2(r* - r) 

dr(po{z,rY = -■ 


( 34 ) 
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To complete the scheme, the maximizing parameter r* has to be calculated from (34) and the result r* = 


r*(z, r, r) has to be inserted to (33). 

We solve this equation for 4>Q{z,r) = ^TrlogA4(2:,r) where M.{z,r) = (z — Xq){z — Xq) + with an 
initial matrix Xq 


^iz,r,T) = ^TrlogA4(z,r,) - ^ 

N r 


^TrA4(z,r*)-i = 


(35) 


Eliminating r* from this set of equation we obtain the effective potential (j){z, r, r) for an arbitrary initial 


matrix Xq. From it we derive the limiting density p{z,t) (14) and the eigenvector correlations 0{z,t) (16) 


p{z,t) = -d:,g(j){z, 0 ,T), 

TT 

0{z, t) = ^ lim {dr(l){z, r, T)f . 

47r r—)-0 


(36) 

(37) 


For the initial condition of the form (f>o{z,r) = ^TrlogA4(z,r) we arrive, after differentiation and some 
algebraic manipulations, at 


r) = - _^ det TrM-^n \ ^^r 

’ '' V7rTrA4-2 -TrA4-V Tr(z - Xo)A4-2 ' + '' 


Ntt' 

0{z,t) = 


.. (a4-1[A4-1;z-Xo](z-4 


(38) 

(39) 


where A4 = A4(z, r*). In the final formulas we set r = 0 to project the results to the z-complex plane. The 
equation for r, (|3^ simplifies for r = 0 to 


^TrA4(z,r.)-i = -. 

N T 


Equations (38), (39) are valid inside the boundary given by 


lTrA4(z,0)-i = i, 


(40) 


(41) 


which corresponds to r» = 0. Outside this boundary 0(z,t) = 0 and p{z,t) = 0. For a normal initial 
matrix Xq, the second term in the spectral density (|^ drops out since [A4“^,z — Xq] = 0 . 


By inspecting the boundary equation (41) one finds a surprising connection to the so-called pseudospec¬ 
trum [31j . a mathematical concept generalizing the notion of the eigenvalue spectrum. Pseudospectrum of 
a matrix A is defined as an subset cr^i of the complex plane z such that 


:-A)- 


1 

> 

€ 


(42) 


where the symbol 11 • 11 is some arbitrary matrix norm and e is the parameter of the pseudospectrum. In the 
e —>■ 0 limit, one recovers the standard eigenvalues as poles of the resolvent (z — A)~^. 

For the initial matrix A = Xq, the boundary of the pseudospectrum subset axg is exactly the eigenvalue 
boundary (41) with ^ and a Frobenius norm ||X||_f = VTizW. From this simple identification 

we conclude that a diffusion model with an initial matrix Xq is also a probabilistic realization of the 
pseudospectrum for the same matrix. 

We finish this section by discussing an equation for the gradient v = in the inviscid limit. The 
equation is equivalent to the one for the potential (j) that we discussed above but in some situations the 


equation for the gradient is more handy to use. The inviscid version of the (28) reads 


d-rV = d,. 


( 43 ) 















It is an inviscid Burgers equation in 2 + 1 dimensions. A general solution to this equation for smooth 


(differentiable) initial conditions can be deduced from the Hopf-Lax formula (24) for the effective potential 
(j) for iV —)■ 00 . The maximum in (24) is achieved for w' = ic* fulfilling the standard extremum condition 


0 = vo{z,w^) - 


w* — w 


(44) 


which is equivalent to 


w = Wt, — TVo{z, w*), 


(45) 


for which v{z,w,t) = vo{z,Wi,). Inserting ru* (45) to this equation we get a solution to (43) 


= vo(z,w + tv), 


(46) 


which is given by an implicit equation for v = v{z, w, r) depending only on the initial condition vq{z, w). The 
parameter w* can be viewed as a labeling parameter for the family of the characteristic lines. These lines 
start to cross when the labeling fails to be bijective i.e. = 0. This singularity condition defines a caustic 
surface - a boundary across which an ambiguity of the solutions arises. The development of multivalued 
solutions is an unwanted feature of the inviscid Burgers equation and is circumvented by constructing shock 
lines along which one cuts the characteristics rendering the solution unique. 


The two-dimensional Burger’s evolution (43) can be simplified in our case due to the rotational symmetry. 
We are interested only in the solutions which depend on the modulus r = |r(;|. In this case the vector velocity 
field is a central field v = with v = |n|, and the vector equation (43) reduces to a scalar equation for 
the modulus of the velocity field 


dT-v = vdrV, (47) 

with a solution given by 

V = VQ{z,r + tv). (48) 

From this solution we can reconstruct the full 2d-solution; v{z,w,t) = '^v{z,r,T) with r = |r(;|. 

6. Examples 

In this section we discuss three examples: (1) the canonical Ginibre evolution for which the initial matrix 
is Xq = 0, (2) the spiric evolution for which Xq = diag(—a,—a...,a,a...) with an equal number of ±a and 
(3) an evolution initiated from a non-normal matrix: (Xo)ij = aSij-i. This matrix has eigenvalues equal 
zero as the initial matrix in the first example but it is not a normal matrix. 

The first example serves as a proof-of-concept. We solve the Burgers equations by the method of charac¬ 
teristics to obtain the spectral density, the eigenvector correlator and the potential function in the large N 
limit. The discussion is accompanied by Appendices B,C where we consider finite w results and a relation 
of the averaged extended characteristic polynomial to the two-point kernel of the underlying determinantal 
process. 

The second example illustrates an evolution of the eigenvalue density initiated from two disconnected 
eigenvalue ’’islands” which grow in the course of time to collide at some critical time. We discuss a novel 
universality arising in the vicinity of the collision. 

The third example demonstrates the dependence of the evolution of the eigenvalue distribution on the 
initial information which goes beyond the eigenvalues themselves. 
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Figure 1; A numerical simulation of the spectral density with an initial matrix Aq = 0 at time slices r = 0.1, t = 0.2 and 
r = 0.5 respectively, an ensemble of 6 matrices of size iV = 1500. Black curves are the large N spectral boundaries and the red 
dot indicates an arbitrary spectator position z where the evolution is probed. 


6.1. Ginibre evolution 

The evolution is initiated from the matrix Xq = 0. The evolution of the eigenvalue density is shown 
in Fig. [2 The spectral density forms a circular eigenvalue ’’island” expanding in time. For each time 
T the ensemble of matrices in this evolution is equiv alen t to a Ginibre Ensemble with a time-rescaled 
dispersion. The initial condition for the determinant is Dq{z, w) = (|zp + I 
potential (18) (j)o{z,w) = log(|zp + Iwp), for the velocity vo{z,w) = w/{\z\^ + |w 


V N 


for the effective 
), and for its modulus 


i^oiz,r) = r/{\zY + where r = |r(;|, respectively. We solve the inviscid 2-1-1 Burgers equation (43) using 
the method of characteristics. Characteristics (45) for this equation are shown in Fig. In the left panel 
we show a plot in 2-|-l dimensions and in the right one its 1-1-1 dimensional section. The evolution takes 




Figure 2: The characteristic lines at z = 1 for the vector 2-1-1 Burgers equation (left) and a 1+1 section (right). The caustic 
cone-like surface on the left plot is denoted by dashed line on the right one. The shock is located on the cones’ axis It start 
from its apex. It is shown a green double line on the vertical axis in the right chart. The vertical line corresponds to r = 0, 
that is to the place where the quaternionic pair (z,in) reduces to (z,0) that lies on the z-complex plane. 

place in the w complex plane however the position on z-plane defines the initial condition. The z-variable 
acts as a spectator on the eigenvalue plane ’’observing” the evolution in the perpendicular w-direction. We 
identify the cone-like caustic surface (left plot) whose apex is located at r = 0 and at critical time Tc = \z\'^. 
This surface is the boundary along which the characteristics start to cross, making the Burgers solution 
multi-valued. This ambiguity develops for r > Tc as the expanding eigenvalue boundary ’’swallows” the 
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spectator at zg- We find the position of the shock line as a locus distanced equally from the caustic surface 
at each given time (the Rankine-Hugoniot condition). Since our problem is radially symmetric, the shock 
is positioned exactly at r = 0 and starts from the critical time Tc- Therefore although the dynamics takes 
place in the whole r space, the shock is always confined to the ’’physical” r = 0 region. Moreover if we 
confine it to r = 0, the spectator at z stays on the shock line at every time t > Tc- We conclude that inside 
the bulk of the spectrum (i.e. the non-holomorphic sector), the observer is constantly on the shock line. 

Having identified the positions of the shock we can write down a solution of (48) for the reduced 1+1 
Burgers equation for our initial conditions 1 ^ 0 ( 2 , r) = r/{\z\^ +r^). It reads 


r -\- Tu 


(r + Tv)^ + ' 


(49) 


It is an implicit algebraic equation for v = v{z,w,t). It can be rewritten as a cubic equation. Solutions 
for different r is plotted in Fig. Rather than showing the solution for the modulus v we show a cross 
section of the vector field which has two symmetric branches ±z/. The right plot shows the solution inside 
the caustic surface with an unphysical branch depicted as dashed lines. For r = 0, that is on the z-plane, 
we have 


v[z,t) = 


0 




for 

for 


r < |z( 
r> kl' 


(50) 


The boundary |zp = t is found by the sewing condition of zero and non-zero solutions. From v we readily 



Figure 3: The solution |ut| for fixed jiq = 1 a-rid three different times r = 0.5, r = 1 and t = 1.5 respectively. It shows the 
development of nonzero solution for t > and the emergence of an additional unphysical solution depicted by dashed lines. 


obtain the eigenvector correlation function (|15|) in the large N limit 

1 


0{z, t) = —v^ = 


0 for T < |z|2 

W(r-k|2) for r>|z|2 


(51) 


We can also find the second gradient g = dz4> of the effective potential. To this end we can use the equation 
(29) which in the limit N —>■ 00 simplifies to drg = dw\v\'^. For our initial condition gQ{z,r) = z/(|zp + r^) 


it gives 


giz,T) = 


Xjz for r < 
zjr for r > 


Now we can use the Gauss law (10) to obtain the limiting eigenvalue density 

p{z,t) = —6 »(Vr- \z\). 


(52) 


(53) 


It is given by a uniform distribution on a disc of radius i/t. For r = 1, results ( 5l|53 ) reproduce results 
of [MlUS], respectively. We could alternatively find the same formulas by directly applying (38) and (39). 


11 





















The maximizer r* can be found by solving the constraint (40): 


0 

\/^ 


We find 




ln|z 

N 


for 

for 


T < |Z| 

r > fol 


for T < \z\ 


Inr + -1 for t > \z\ 


( 54 ) 


(55) 


in which as the effective potential (j) and its gradients are 


The last formula was obtained for t = 1 in [5^ . 

So far we have discussed the limit N ^ oo, 
equal to the electrostatic potential and the quaternionic Green’s function (30). What about finite N7 For 
finite N the order of calculating the average and the logarithm in 0 and (|l8[) matters , so ()) is at best 
only an approximation of $ for large but finite N. So clearly our method is not able to get insight into the 
finite N corrections. Surprisingly as we show below, the behavior of the diffusion kernel D{z,w,t) on the z 
complex plane, that is for w = 0 reveals the same type of finite size effects as known from exact calculations 
of the eigenvalue density for Ginibre ensemble [24] . For our initial conditions the diffusion kernel (25) is 




(56) 


On the z-complex plane that is for r = 0 this integral significantly simplifies and it can be for large TV 
calculated using the saddle point method. There are three saddle points 

r(, = 0,r± = \z\^. (57) 

When r approaches |zp from above, the two points r'j. approach each other along the real axis to coalesce 
for r = |zp at r' = 0. When r becomes smaller than |zp the two points move on the imaginary axis. Near 
the critical value r = |zp it is convenient to introduce a rescaled parameters 

r' = dAr-i/4, 1^1 = ^ + ^^-1/2 (58) 


and use them in the calculations. One finds that for large N the diffusion kernel behaves as 


D{z = y/T + r]N ^/^,r = 0,r) 



(59) 


where erfc is the complementary error function. This result holds not only on the z-complex plane but also 
sufficiently close to the z-complex plane that is for r approaching zero as or faster: r = 

Indeed in this case the argument of the function Ig is a finite number and this function behaves as a 
constant for large N. The error function behavior has the same form as the finite size expression for the 
eigenvalue density function. This indicates that there may be some deeper relation between the diffusion 
kernel (averaged extended characteristic polynomial) and the two-point kernel known from the considerations 
of the underlying determinantal process for Ginibre matrices. We work out this analogy at heuristic level 
in Appendices B,C. 


6.2. Spiric case 

We now consider diffusion initiated from a diagonal matrix Xq = diag(a,..., a, —a,..., —a) with the same 
number of a and —a eigenvalues. In the course of time two initial eigenvalue ’’islands”, initially concentrated 
around ±a, expand to collide at some critical time as shown in Fig. [^ The initial condition corresponds 
to 

Do{z,r) = (r^ + \z - ap)^ (r^ -f |z -f a\‘^)^, 
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Figure 4: The spectral density dynamics for an initial matrix Xq = diag(l,1, —1,—1) before, at and after the critical 
time Tc = 1, the ensemble consisted of 6 matrices of size N = 1500. The green, red and black dots denote three observers 
useful in the analysis of the evolution. 



Figure 5: Characteristic lines in the spiric case for fixed observers zq = 0.5 (black dot on Fig. [^, zq = 0 ((green dot on Fig. 
and Zq = 1 i (red dot on Fig. with a = 1. From left to right, the shock formation occurs later in time which corresponds 
to meetings between observers zq and the expanding spectral boundary. 


or if we write it for the modulus of the velocity field (471 

1 


Mz,r)= ^ 


2 + |z — a\ 


(60) 


In Fig. [^we show charecteristics for three different values of the observer position z, including z = 0 where 
the collision takes place. As we can see the dynamics in the r plane behaves qualitatively in the same way as 
for the Ginibre case. In particular , characteristics form the same type of caustic surfaces with a shock line 
present after some critical time Tc = yypqriyp- This critical time can be found by sewing the solutions or 
by finding the position of the caustic surface apex as a function of z and t. In Fig. [^we show a comparison 
of the caustic surfaces for Ginibre and spiric evolution for r = 0.5. For this inititial condition the solution 


(481 of the inviscid Burgers equation takes the form 


2v = 


r + TV 


r + TV 


z + a\^ + {r + Tv)^ |z — op + (r + 


(61) 


For r = 0 the solution can be written explicitly 


v{z, r = 0;a) 


0 for z ^ 5 

y \/^ - l“P “ + 5‘5'a for z €S 
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Ginibre 


Spiric 




Figure 6: Caustic surface in the variables (r, z) for Ginibre case (left plot) and spiric case (right plot) with time fixed at 
TO = 0.5. 


where Sa = \/t‘^ + 4^^, Za = za + za. The symbol S stands for the interior of the spiric section defined by 
the contour 


r(|a|2 + |zP) = 


(63) 


Generally spiric section is a curve obtained by intersection of a torus and a plane parallel to its rotational 
symmetry axis. We plotted spiric curves as contours around the scatter plots for eigenvalue densities in Fig. 
1^ The eigenvector correlation function inside the spiric section S is 

.,2 


0{z,t) = — = 


N ( T 


TT r^TT V 2 


1 


- a - Z 


S„ 


(64) 


We can also calculate the diagonal element of the Green’s function using the equation (29) that in the limit 
fV —)• oo simplifies to drg = dzV^. For the initial condition as 50 = 5 + r^+|z-ap ) solution 

reads 


g{z,r = 0,t) = 


2 _ aSq 

T 2TZa 


for z ^ 5 

+ c{z, a) for z £ S ’ 


(65) 


with a constant c{z, a) = obtained by the sewing condition along the spiric section S. We use the Gauss 
law to obtain the spectral density as 


p(z,t) = 


2|„|2 


Sa{2 - r|aF) + T^\a 

2TTTZ^Sa 


( 66 ) 


The same results can be obtained from the calculations of the effective potential (j) using the equations (38), 

f 0 for z ^ 5 


(391 and (40) for r* 




- lap - Izp + iS'a for z e S 


Outside S we get 


and inside S 


h°^\z,T) = iln|z - ap + ^ln|z + ap, 


‘(^,r) = iln(^ + I^.)+l(|«P + |zP-I-l5. 


(67) 

( 68 ) 

(69) 


14 













Pmarg(Xi^—*1) 


Omarg(^-^— 




Figure 7: Numerical results in the spiric case for a = 1 and critical time r = 1 for eigenvector correlator (right plot) and 
spectral density (left plot) averaged over the imaginary axis, the ensemble consisted of 6 • 10^ matrices of size N = 1000. 


We recall that Sa = \/r^~T4^, Za = za + za, as defined after ( [d^ . In Fig. [^we show a comparison ol' 
theoretical predictions for the limiting eigenvalue density (661 and the eigenvector correlation function (641 


with Monte-Carlo simulations. The agreement is very good. 

We complete the discussion of the spiric case by deriving a finite N formula for the diffusion kernel 
D{z,r,T) 


D{^,r.r)= f r 'exp (-iV^yf) 


h ) (/2 + |2 _ o|=)T(r" + |e + ap)Td/. (70) 


We are interested in the behavior close to the origin x = 0, r = 0 for r close to the collision time. Without 
loss of generality we perform calculations for a = 1. In this case the collision time is r = 1. We set r = 0 
and zoom into the vicinity of the critical region 


r' = z = 


T = 1 + tN-^/^ 


(71) 


where the saddle points merge. After expanding the solution (70) we obtain an asymptotic formula 
D{z = r = 0,T=l + r 


128V 




- iv + vf -t) 


(72) 


which is plotted in Fig. Again, the result holds not only for r = 0 but more generally for r = 0{N 
6.3. Non-normal Ginibre case 

We consider now diffusion initiated from a matrix {Xo)ij = aSij-i. The matrix Xq has all eigenvalues 
equal to 0. The initial eigenvalue distribution coincides with the one for the Ginibre case but afterwards 
the eigenvalue density obeys completely different evolution. Three snapshots of this evolution are shown in 
Fig. § The initial distribution concentrated initially at zero instantaneously expands to a circle of radius 
|a|, which corresponds to the pseudospectrum of the matrix. In the course of evolution the support of the 
density takes the form of the growing annulus. After a finite time r = |ap the inner radius of the annulus 
shrinks to zero and eigenvalues fill up a full disk. 


Let us show this by direct calculations. The matrix A4 = (z — Xo)(z — Vq)^ -\-r^ (351 has for our choice 
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Figure 8: Evolution of the characteristic polynomial in the vicinity of the collision for rescaled times t = —l,t = 0,t = 1. Black 
contour on top is the large N boundary of non-zero spectral density obtained from | |63[ l. 

2 
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0 

-1 

-2 

- 2-10 1 2 - 2-10 1 2 - 2-10 1 2 

Figure 9: Evolution of the spectral density for non-normal initial condition (Xo)ij = with time snapshots at t = 0.2, 

r = 0.5 and t = 1, the ensemble consisted of 6 matrices of size N = 1500. 





of Xq a tridiagonal form 


/a 

b 



b 

a 

b 


0 

b 


b 

0 


.. 0 \ 

0 

b : 
a b 
b d) 


(73) 


with a = Izp -I- -I- |ap, b = —za and d = |zp -|- r^. The determinant of this matrix can be calculated 

explicitly: det At = ^ (d(a+ — a.^) — where A = — 4|6p and a± = | (a ± A). The 

initial effective potential for large N is 

(l)o{z,r) = ^TrlogAt « Ina+ = In ^ -I- a/o^ - 4|zp|Q;p^ (74) 

where a = \z\'^ + + \a\'^. We have neglected 1/A terms which disappear in the limit A —t oo. The value 

r* which maximizes the expression in the Hopf-Lax formula (26) can be calculated from the equation (34). 
For r = 0 we get 


r 0 for z ^ A 

( for z £ A 


( 75 ) 
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where we used the notation Tq, = \/t^ + 4|ap|zp. The boundary of the annulus A is given by the radii 


l^li = \/|aP ± T. 


(76) 


We see that at the beginning of the evolution, r « 0, the annulus is infinitely narrow forming a one¬ 
dimensional circular pseudospectrum. On the other hand for r —> |ap the inner radius tends to zero and 
the annulus becomes a disk. Inserting r* to (26) we obtain the potential 

|2 


In |a| 
In UP 


for 

U| e (0;U|_) 

- ^ for 

T 

U| e (U|_; \z\+) 

for 

U| e (U|+;oo) 

is obtained 

from the Poisson 


(77) 


differentiation the effective potential twice 




(78) 


on the support of the annulus and zero otherwise. Using (39) we obtain the eigenvector correlation function 
to the support of the annulus 


1 


0(ur) = ^(r„-Up-|ap). 


(79) 


In Fig. we compared the prediction of the two formulas given above with numerical simulations. The 
agreement is good. Deviations are observed only close to the boundaries and can be attributed to finite-size 
effects. It is easy to check that all the expressions transform to the corresponding expressions for the Ginibre 



Figure 10: Numerical analysis of the eigenvector correlator (right plot) and the spectral density (left plot) for the non-normal 
initial condition, the ensemble consisted of 3 ■ 10^ matrices of size N = 1000 with parameters a = 1 and time t = 0.8. 

evolution when ct —> 0. 

Let us shortly discuss finite N effects for the averaged extended characteristic polynomial. We are 
looking for an universal behaviour near the origin for t — i 0. Around this space-time point an instantaneous 
transition from N eigenvalues positioned at the origin to the ring of radius |a| happens. In our example the 
characteristic polynomial is given explicitly as 

2N /2Nrr'\ 

D{z,r,T) = - J r'expf—° - j det M{z,r') dr' (80) 
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We consider the following scaling around the origin 

r' = |z| = 

|a| = cciV-i/®, T = 


(81) 

(82) 


Inserting these formulas to the equation above we find an asymptotic function for large N and for |a| = 




(83) 


We see that it reveals a non-perturbative divergent character near x = 0. For the normal Ginibre case we do 
not find the second term and therefore we identify it as a consequence of the non-normal initial condition. 
Concluding this section, we note that the formula (781 was obtained in |22j for an initial matrix 

g277j(JV-l)/JV^^ 


= |a|diag(l,e' 


27T’i/N AnilN 


(84) 


with N initial equidistant eigenvalues lying on a circle of radius |a|. This result is not surprising since the 
matrix Xq is unitarily equivalent to a circulant matrix 


(-^0 )^i ^ (^o)ij 4“ 1; 


(85) 


which in turn differs from the non-normal initial matrix considered in this paper by just one element whose 
effect can be neglected in the large N limit. We also mention that finite N spectral densities and large 
deviations have been recently studied in [33]. 


7. Spectral stochastic equations 

An interesting issue arises, if the stochastic dynamics of the elements of non-hermitian random matrices 
corresponds to some stochastic equation for the spectra of the underlying matrix. Such phenomenon holds 
in the case of Gaussian Unitary Ensembles (GUE), as well as in the case of Circular Unitary Ensembles 
(CUE), as shown in the seminal paper by Dyson |5]. In the simpler case of GUE, the corresponding Langevin 
equation for real eigenvalues reads 


dt 


dAi — dBi + -- 




Xi X-i 


( 86 ) 


where dBi reflects the Brownian dynamics and the second, ’’drift” term comes from the Van der Monde 


determinant. We neglect the optional Ornstein-Uhlenbeck term —Xt dt on the r.h.s. of (86) since it only 


contributes to freezing the diffusing end-points of the spectra in the stationary limit. The corresponding 
Smoluchowski-Fokker-Planck equation written for the resolvent G{z,t) of the spectral density p{X,t), takes 
(after rescaling the time r = Nt and in the large N limit) the form of complex, inviscid Burgers equation, 

i.e. the eq. ([^ piM]. 

Equations for non-hermitian Gaussian ensemble obtained in this work exhibit structural similarity to the 
hermitian ’’Burgulence”, hence a question arises, if some stochastic dynamics for complex eigenvalues does 
exist as well. Recently, Osada m had examined a stochastic two-dimensional Coulomb system defined as 


dA, = dB, 


( 2 ) 


E 


dt 

Xi Xi 


= dB. 


( 2 ) 


-V dt 


(87) 


(o') 

where di?) is a two-dimensional (complex) Brownian walk. In particular, he has shown that such interacting 
Brownian motion equipped with trivial initial condition (all Xi put to zero at t = 0) leads to limiting 
distribution of Xi representing a uniform disk, therefore resembling the Ginibre Ensemble spectrum. 
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From the point of view of our analysis this result is curious, since we have proven that the dynamics 
of eigenvalues is intimately connected to the dynamics of eigenvectors, which seem to decouple completely 
from “hypothetical” eigenvalues of Ginibre ensemble in (87). We suggest the resolution of this puzzle. 

It is well known, that in so-called Normal Random Matrix model [35] with axially symmetric potentials all 
correlation functions can be expressed in terms of holomorphic functions of a single variable. Moreover, the 
exact integrability of such models can be linked to (2-|-l)-dimensional Burgers equations. Quite remarkably, 
in the case of the potential V(z,z) = |zp, the correlations of the Normal Random Matrix model are identical 
to the correlations of the Ginibre Ensemble |36|. We therefore conjecture, that the stochastic equation (871 
corresponds rather to Gaussian Normal Random Matrix model than to the Ginibre Ensemble. 

Because normal matrices are diagonalizable by a single unitary transformation, alike the hermitian 
matrices, the eigenvectors decouple from the eigenvalues which explains the lack of these degrees of freedom 
in (871 and its formal similarity to hermitian stochastic equation (86). Based on this and the significance 
of eigenvectors presented in this paper we find it highly probable that the Ginibre-like dynamics of (87) is 
accidental and proper stochastic equation governing the dynamics of non-hermitian random matrices is not 
known. 

The speculative links between our approach for generic complex matrices and non-Gaussian Normal 
Random Matrix models represent a challenge, which we plan to address in the forthcoming publications. 


8. Conclusions 

We have shown that a consistent description of non-hermitian Gaussian ensemble requires the knowledge 
of the detailed dynamics of co-evolving eigenvalues and eigenvectors. Moreover, the dynamics of eigenvectors 
seems to play a superior role (at least in the N —>■ oo limit) and leads directly to the inference of the 
spectral properties. This is a dramatically different scenario as compared to the standard random matrix 
models, where the statistical properties of eigenvalues are of primary importance, and the properties of 
eigenvectors are basically trivial due to the their decoupling from the spectra. We have considered examples 
of Ginibre, spiric section and non-normal Ginibre where the formulas for spectral density, the 1-point 
eigenvector correlation function, electrostatic potential and universal functions were obtained and positively 
crosschecked with numerical simulations. By studying the dynamics of characteristics we anticipated the 
novel universality obtained for the spiric case as an error function type. We obtained compact formulas for 
both spectral density and eigenvector correlator for which the latter unraveled a promising determinantal 
structure. The diffusion equation, as an equation exact for finite N, was used mainly to obtain the universal 
behavior. 

We conjecture that the hidden dynamics of eigenvectors discovered by us and described for the Gausssian 
non-hermitian ensemble, is a general feature of all non-hermitian random matrix models, and has to appear 
systematically in 1/fV expansion. 

Our formalism could be exploited to expand the area of application of non-Hermitian random matrix 
ensembles within problems of growth, charged droplets in quantum Hall effect and gauge theory/geometry 
relations in string theory beyond the subclass of complex matrices represented by normal matrices. 

One of the challenges is an explanation, why, despite being so different, the Smoluchowski-Fokker-Planck 
equations for hermitian and non-hermitian random matrix models exhibit structural similarity to simple 
models of turbulence, where so-called Burgers equation plays the vital role, establishing the flow of the 
spectral density of eigenvalues in the case of the hermitian or unitary ensembles and the flow of certain 
eigenvector correlator in the case of non-hermitian ensembles. 

We believe that our findings will contribute to understanding of several puzzles of non-hermitian dynam¬ 
ics, alike extreme sensitivity of spectra of non-hermitian systems to perturbations We also hope 

that the quaternion extension used in our paper may help to understand better the mathematical subtleties 
of the measure of non-hermitian operators [37j . 
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Appendix A. Derivation of the diffusion equation (201 


In this appendix we demonstrate that the averaged extended characteristic polynomial (19) obeys the 
diffusion equation (201. The determinant in ( |l^ can be expressed with help of Grassmann variables as 

D{z,w,t) = J V[X]V[ri,^]P{X,T)expTGiX,z,w,r],^), (A.l) 

with the object Tq given by 

Tg{X, z, w; ??, ^) = ^ - iyij {fj.rjj - ^ {zfj.rji + z^i^i + , (A.2) 

where and ^ are Grassmann variables and Xij = Xij + iyij- With the help of heat equation (§ for 

the joint probability density function P{X,t) we obtain 




)-) 


exp Tg = 


drD{z, w,t) = j V[X,y, (\ {drP) exp Tg = X /" T>[A, 77, | ^ {c 


where we integrated by parts twice and we used 

expTc = {fji-rjj + expTg = expTc, 

dl,. expTc = (fjirjj + expTc = expTc. 


On the other hand we have 


d^uw exp Tg = '^ exp Tg, 




from which it follows that 


^d^^D{z,w,T) = ^ 




(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 


We see that the expressions on the right hand side of (A.31 and (A.7) are identical and thus we have 
^dwwD{z,w,T) = drD{z,w,T) (20). As a side remark we note that this calculation can be almost verbatim 
repeated for the averaged extended ’’inverse” characteristic polynomial 

which obeys the diffusion equation in the opposite time direction r —> —t 

1 


(A.8) 


— f) P — _ f) - P 


(A.9) 
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Appendix B. The r ^ 0 regime as a Wishart/chiral deformation 


The majority of results discussed in this paper are confined to the ’’physical” region where r —0. One 
can however keep the parameter r nonzero. In the Coulomb gas interpretation, this deformation introduces 
a complex nonlinear interaction between the eigenvalues of unknown interpretation. 


We consider the Ginibre case. The expression in the Hopf-Lax equation (26) 


(/)(z, r, r) = max ln(|zf + r'"^) — 


(r' - rf 


is maximized by r* obeying a cubic equation 


(r* - r){rl + = r*T, 


(B.l) 


(B.2) 


The corresponding spectral density p{z,r,T) and the eigenvector gradient v{z,r,T) are shown in Fig. B.ll 


There are no critical points for r ^ 0 and instead we see a smooth crossover between two phases. This is 
in accordance with the fact that the shock is present only on the r = 0 plane. Consider the formula for the 



Figure B.ll: The large N spectral density (left plot) and eigenvector correlator \v-t-\ (right plot) for r = 0.05 (blue line) and 
r = 0 (red line). 


averaged extended characteristic polynomial for the Ginibre evolution 

' 2Nrr' 


2N / j 

D{z,r,T) = - / r'exp ( -N- 

Jo V 


(1^ 


dr'. 


(B.3) 


Applying the Newton binomial formula to (|zp +r'^)’^ an using the following integral representation of 
Laguerre polynomials 


dqq 




\ 




N 


-r I e ' 


we get 


D{z,r,T)=e ( 

k=0 ■ ^ 


Nr^ 


N\2 


(B.4) 


(B.5) 


It is not surprising that Laguerre polynomials show up in this context since the the Ginibre ensemble is 
closely related to the Wishart/chiral ensemble where they naturally occur. In fact by recalling the definition 

pi 


D{z, r, r) = (r^ + {z - X){z - X^))^ 


(B.6) 
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we see that for z —)■ 0 the object within the brackets is a Wishart matrix and 


D{z = 0, r, r) = 


For r —>• 0 we have D{z = 0, r = 0, r) = A^! (^) 






(B.7) 


N 


Appendix C. The kernel structure 

We argue that in the case of Ginibre matrices, the characteristic polynomial D in the r —0 limit 
has essentially the same information as the microscopic kernel of the underlying determinantal process. 
This stems from the observation made in |39| where the author connects the n-point correlation function 
of a Ginibre matrix model to a random matrix QCD partition function with n quarks with appropriately 
decreased matrix size: 


det(AAr(zj,Zj))ij=i..,„ = c(zi, ( ]^det(zi - A)det(zj - A'^') 


(G.l) 




where det(AAr(zi, ..n = (nt Tri5^(zi — is the correlation function averaged over N dimen¬ 

sional matrix and c denotes a known z dependent proportionality factor. 

We study closer the n = 1 case because for this parameter the prefactor c is the weight function w{z) 
and the averaged determinants are exactly the characteristic polynomial D in the r —> 0 limit. Therefore, 
in this particular case, we obtain a on-diagonal kernel formula 


Kn{z,z) = Cnw(z)D^^ ^\z,r = 0 ,t), 


(G.2) 


with some numerical constant Gjv. The off-diagonal kernel is not so easily obtainable by the above formula 
however we make an educated guess based on the symmetry of arguments and the result of Akemann and 
Vernizzi 135]. We write the full kernel as 


Kn{z,v) = Cn^/w{z)^/w{v)D^^ ^)([z,t>],r = 0,r), 
where the two-argument characteristic polynomial D was created by substituting |zp 


zv: 


r.(N-i)n , , 2N f ,y+r^\rf‘^Nqr\,^^ 

D'' >{[z,v],r,T) = — J qexp\^-N --- + zv) 




dq, 


which for r = 0 it is readily solved as 


^ ^ fNzvY 1 
W 


(G.3) 


(G.4) 


(G.5) 


fe =0 


To complete the argument demonstrating that the full information resides in the characteristic polynomial 


D we will find an a priori unknown weight function w{z) present in formula (C.3) from the D alone. This 
is done by using decomposition in the biorthogonal basis 


N-l 

^ CmXu{z)Yi{v), 
k,l—0 

where Xi,Yi are biorthogonal polynomials with respect to an unknown measure W 


(fzW[z)Xi{z)Yj{z) = gij. 


(G.6) 


(G.7) 
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The expansion terms Cij and the bi-orthogonality matrix gij are inverses gij = c~^. We proceed by finding 
a basis Xi, Yj in which Cki are diagonal. Then we infer the formula for W by considering its moments (C.7). 
We start with a formula (C.5) which is already in a diagonal form with 


Cfcfe - Ur 


N-l-k 


Xfe(z)=z^ Yk{v) = v^ 

The moments are therefore given by 

/= I; (7) 

By assuming the radial symmetry and setting |zp = p we obtain 

^0 

from which the characteristic function is given by 

POO 

Mw{t) = / dpW{,/^)e^*P = 

Jo 


N\ / T \ N-l-k 


! \n) 


ttNI (iV/T)fc+i ■ 


(N/t)^ ^ itr 


TriV! 


N 


-1 


It is exactly the characteristic function of an exponential distribution Ae with A = N/i 


(C.8) 


(C.9) 


(C.IO) 


(C.ll) 


(C.12) 


This procedure gives the unknown weight w{z) = e 
1 ( 

tzN\ \ T ) 


along with the normalization coefficient Cm = 
W{\z\^)=Cnw{z). (C.13) 


The kernel is therefore equal to 


1 / N \ ^ 

Kn{,z,v) = —ei^pi- — {\z\'^ + \v\'^)\ Y, 


f Nzv\^ 1 


\ 


k\' 


(C.14) 
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